A Modified Mualem–van Genuchten Formulation for Improved Description of the Hydraulic Conductivity Near Saturation

The unsaturated soil hydraulic properties are often described using Mualem–van Genuchten (MVG) type analytical functions. Recent studies suggest several shortcomings of these functions near saturation, notably the lack of second-order continuity of the soil water retention function at saturation and the inability of the hydraulic conductivity function to account for macroporosity. We present a modified MVG formulation that improves the description of the hydraulic conductivity near saturation. The modified model introduces a small but constant air-entry pressure (hs) into the water retention curve. Analysis of the UNSODA soil hydraulic database revealed an optimal value of24 cm for hs, more or less independent of soil texture. The modified model uses a pressure dependent piece-wise linear correction to ensure that deviations between measured and fitted conductivities between pressure heads of 0 and 240 cm were eliminated. A small correction was found necessary between 24 and 240 cm, and a much larger correction was needed between 0 and24 cm. An average RMSE in logKof only 0.26 remained for a data set of 235 samples. The resulting modified MVG model was found to have small systematic errors across the entire pressure range. The modified model appears well suited for large-scale vadose zone flow and transport simulations, including inverse modeling studies. MANY VADOSE ZONE flow and transport studies require description of unsaturated soil hydraulic properties over a wide range of pressure heads, h. The hydraulic properties are often described using the poresize distribution model of Mualem (1976) for the hydraulic conductivity in combination with a water retention function introduced by van Genuchten (1980). The soil water retention equation, u(h), is given by u(h) 5 5 ur 1 us 2 ur [11 jahj] h # 0